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ABSTRACT 

A convenient procedure Is described to determine the thermal behavior 
(thermal expansion coefficients and thermal stresses) of anpleplled fiber com- 
posites using a pocket calculator. The procedurt* consists of equations and 
appropriate graphs for various (-H)) ply combinations. These graphs present re- 
duced stiffness and thermal expansion coefficients as functions of 40 In or- 
der to simplify and expedite the use of the equations. The procedure is appli- 
cable to all types of balanced, symmetric fiber composites Including Interply 
and Intraply hybrids. The versatility and generality of the procedure Is Il- 
lustrated using several !?tt*p-by-step numerical examples. 


l.t) INTRODUCTION 

Themal expansion coefficients and thermal and residual strains and 
stresses In angi •plied laminates (fig;;. 1 and 2) are frequently required for 
the Initial sizing of structural components made from flbt'r composites. These 
coefficients strains end stresses ;ire commonly referred to as composite thermal 
behavior. The significance of composite thermal behavior particularly lamina- 
tion residual stresser, are extensively discussed in reference ]. Thermal ex- 
pansion coefficients, thermal and lamination residual strains and stresses are 
determined using composite mechanics and laimlnate theory usu;ill.y available In a 
computer code (refs.. 2 and 3). A computer code was used effectively (ref. 4) 
to evaluate lamination residual stresses In angleplled laminates- and thereby 
assess the effects of these stresses on the structural intei’grlty of composites. 
It is generally recognized that the use of a computer code Is expedient and 
quite general. However, it does not provide the user with Insight and instant 
feedback of the laminate thenruil behavior and capability as he proceeds with 
the design/analysis of the component. Also, a computer code may not be readily 
accessible to tlie user. 

A convenient procedure (method) is described in this paper which can be 
used to determine tbe thennal behavior of angleplled laminates. The procedure 
is suitable for hand calculations using a pocket calculator. Tt consists of 
simple equations and appropriate graphs of (40) ply combinations from the most 
frequently used composites (figs. 3 to 18). Graphs for other composites can be 
gener. ed by using uniaxial properties from table I .and well-knowia transforma- 
tion equations (ref. 4). The procedure makes use of the well-known transforma- 
tion equations, and laminate t ico-y equations. Its structure is similar to 
that in reference 5. The procedure can handle all types of composites includ- 
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Ing InCerply and Intraply hybrids , Th* procedure it llluaCrated ualng selected atep* 
bya tep numerics 1 examples . The ^iicuaa Ion in this paper la 1 1ml ted to linear mechan* 
leal and thermal behavior of cumpoaltes. Several additional examples and Implies* 
tlons to dea Ign app! Icatlona are described In a NASA TP report under preparation. 

The paper nnd numerical exumplea are written mainly in a tutorial manner 
In order to Illustrate the step-by-step procedure. For this reason, the vari- 
ous sections are self contained as much as Is practical at the expense of some 
duplication. The notation used Is defined when It first appears and also sum- 
marized under symbols for convenient reference. Customary units are used 
throughout the numerical examples *lncc these serve primarily to Illustrate 
step-by-step nut.ierlcal calculatlo.is. Appropriate conversion factors are given 
In the symbols. 


2.0 THKRMAI. EXPANSION COEFFICIENTS 


^e ln*plane thermal expansion coefficients (TECs) o^xx *^yy 

sngleplled laminates are determined from the follow! ig equations 




exx 


1 

f 

^exx 


»cyy E 


cyy 


{vpejjQoil " Vcxy‘?021>«Oll (^012 - VcxyQ022)«a2i] 

^ VpoljfQiii - VcxyQt2l)"Ul <*^112 ■ ''cxy‘U22>“t22]} 

{vpfl[(Qo 2 l - VcyxQ 01 l)a 011 + (Qe22 " ycyxQe2 1 )« 022 ] 

''pof(‘U2l ■ Vcyx‘}m)«Ul <‘^122 " ''cyx‘^t21>“t22j| 


( 2 . 1 ) 


( 2 . 2 ) 


The composite 
V and V 


inuuu 1 1 


cxy 


cyx 


'CXX 

are given by 


and the composite Poisson's ratios 


-exx " ^exx 


Qexv. „ ,, Qexy 

Ecyy Qcyy ' 


<cxx 


Qpyy 


vexy ■ t: ' - ' ■* ; Vcyx 


. ^exy 


‘exx 


J 


The "reduced stiffness" Q^'s are given by: 


^exx “ ''p0%ii VpoQtii 

‘^cyy " '^p0^022 '^p0‘U22 

^exy “ '^ po ‘^012 ^ % o ‘^!.12 “ ^cyx 


(2.3) 


(2.4) 


w 


Th« Q 0 *a are obtained froa eh* appropriate figures 3 to 18. The Q£ ■ are 
equal to Qq'b at 6*0 in figurea 3 to 16. Note Q£21 ” <> 6 '* 

and oi^'* *<’• obtained from the aame figurea. The parameter v n denotea thi 

thickneaa ratio of the ^-pliea to the total laminate thickneae vnile V q ia 
the corresponding ratio for the 0 ~pliea. VpQ and VpQ aatiafy the identity 


Vpe ^ Vpo - 1 


(2.5) 


The following procedure ia convenient to calculate numerical valuea 

for a and a for a given [0/4’0l_ APL using equations (2.1) and 
exx cyy *• s 

( 2 . 2 ): 

1. Obtain values for Qon, Qe22» ^812 * 'ie21t Qtll» Ql22» Qll2 " *^121 

O 0 ll» O 02 ?« otll 0122 from the appropriate figures note that ten values 

are needed. 

2. Calculate values for VpQ and VpQ, respectively. 


^ thickness of ■fO-plies 
P® thickness of API. 


( 2 . 6 ) 


thi ckness of 0 -plle s 
thickness of APL 


(2.7) 


3. Qilculnte valuea for Qexx* Qcyy and Qexy using equations (2.4), 
using the information obtained In Item ( 1 ) above and that obtained from equa- 
tions (2.6) and (2.7). 

4. Calculate values for Ep^x* *'*cyy» '^exy ''cyx using the informa- 

tion obtained in Item (3) above. We use the well known relationship 


>'cxy E.. 


( 2 . 8 ) 


to check our numerical values. In suiamary, we need to look-up ten quantities 
from the graphs and calculate a minimum of seven others in order to determine 


the thermal expansion coefficients (TECs) 
( 2 . 1 ) and ( 2 . 2 ). 


and a, 


using equations 


It la worth noting in equations (2.1) and (2.2) that the AIM, TECs depend 
on the properties of the 0± -plies (Qq's), the O-plles (Q^'s) and the APL (in- 
tegrated) properties F,,. and Vq. Also, the APL Poisson's ratios (Vj.) restrain 
the TECs of the APL since these quantities are prectuled by a minus sign and, 
therefore, subtract from the total. Furthermore, the shear moduli contribute 
to the TECa through the Q 0 's. In addition, equations (2.1) and (2.2) can be 
easily extended for more than one set of ±0-ply combinations. Similar terms 
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'^p02l- ^ ^ etc,, lire added to accommodate this case aa will be de- 
acrlbcd with a maarical example later 

Example 2.1. - Calcul,}Ce the TECe of the 8-ply [±30/02] APL made from 

AS/E composite: (1) Following the procedure outlined in Item (1) above we 

obtain the Q'a from figure 7 and the a* a from figure 8 both at 8 ■ ±30 
and 0 (within curve reading accuracy): 



‘^021 • ^012 " 

«011 • 
. 0022 ■ 
(2) Pulltntflng Itei 


11.3x10^ pal 

2.9x10^ pal 

3.7x10** pal 

-2.6x10“** ln/lii/‘*F 
12.8x10“** ln/ln/“F 
1 (2), we calculate 



*^t21 “ *itl2 • 

«tll ■ 
‘>122 ■ 

Vp 9 and Vpo 


18.7x10** pal 

2.0x10** psl 

0.6x10** pal 

0.4x10"** ln/ln/°F 
16.4x10“** ln/ln/°F 


Vp0 • ^/« “0.5 


Vpo “ A/8 - 0.5 

(3) Following Itei.. (3), we calculate the Qc'a using equations (2.4) and 
carrying the unite, aelectlvely for convenience 


*^cxx “ '^p0*J0U '^pO^^eil 

- 0.5 K 11.3x10** + 0.5 X 18.7x10** - 15.0x10** pal 

Qcyy “ VpeQe22 VpoQt22 

- 0.5 X 2.9x10** + 0.5 X 2.0x10** - 2.45x10** pal 


^cyx ■ *?cxy “ 


'^p0‘^ei2 ''p0‘^U2 

0.5 X 3.7x10** + 0.5 X 0.6x10** - 2.15x10^ psl 
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(A) KullowiitR item (4), wv calculatu *''cyy '*c 


:xy 


6.2 


E - Q - . 15.0x10** - . 13.1x10** psl 

cxx cxx 'Vcyy !».A^vin6 


2.45x10' 


^cyy " *Jcyy " * 2.45x10^ - ^ - 2.14x10^ pal 


cxx 


15.0x10** 


. Ssjii 


2.15xlO*» 


Qcyy 2.45x10*’ 


0.878 


. !k2U* 


'cyx q 


cxx 


2. 15x10^ 

15.0x10* 


- 0.143 


Check; '’cyx " '*cxy g 


cyy 


cxx 


2.14x10** 

0.143 - 0.878 r >0.143 - 0.143 O.K. 

IJ.lxlO* 

(5) UHlng the above information In equations (2.1) and (2.2) and cancelling 
10** with 10”* within the braces: 


cxx E 


■XX {'^l’®r**^011 ‘ '’cxy‘^021>“011 ■*■ ^*^012 ” '*cxy^022>“022] 

'^poQ^^ill " '*cxy‘^£21^“Ul *^£12 " '’cxy*^£22)“£22] } 

6 (o.5r(11.3 - 0.878 x 3.7)(-2.6)+(3.7 - C.878 X 2.9)(12.8)1 

13.1x10 I •- -> 

+ 0.5[(18.7 - 0.878 x 0. 6) (0. 4) + (0. b - 0. 878 x 2 .0) ( 16. 4)’]| 

- 7 (o.5r-20.9 + 14.81 + 0.5f7.27 - 19.0]} 

13.1x10** ^ J L JJ 


*cxx 


-0.68x10”** ln/ln/"F 


(Cont'd.) 


“eyy * \''p0l^‘^e21 " ''cyx‘^01l>“0ll * ^%22 ” ''cyx‘J021>“O22j 


^ ■ '’cyx‘UlP“llI ^ <<Jl22 “ ''cvx‘^i2p“«22]} 

^ — X (o.5[<3.7 - 0.143 X U.3)(-2.6)+ (2.9 - 0.143 x 3.7)(12.8)1 

2.14x10° I 

+ O.iQo.b - 0.143 X 18.7U0.4) + (2.0 - 0.143 x 0.6)(16.4i]J 
2 14 ' ' 10 ^ {o.5lj-5.42 + 30 . 3 ) + 0.3Q-0.830 4 31.4]| 


cicyy ■ 13.0x10^ In/ln/^F 


Examp le 2.2. - Calculate the TECa of the [0/t45/90jg MM. made from AS/E 
composite. This APL is cixnmonly called pseudo-isotropic or quasi-isotropic 
laminate meaning that the laminate behaves like an isotropic material with 
respect to its In-plane elastic and thermal properties. The meaning of this 
terminology will also bt> illustrated by a numerical example later, /gain we 
follow the procedure outlined in items ( 1 ) through ( 5 ) and extend it to three 
different ply configurations (0, ±45, 90); 

(1) Obtain from figures 7 and 8 (wit)iin curve-reading accuracy) and using the 
first subscript 0 to denote all three conditions: 


Q and Q 

0 - 445 

0-0 

0 

■ 

Qyjj (IO 6 psi) 

b.l 

18.7 

2.0 

Qe 22 psi) 

6.1 

2.0 

18.7 

^621 • Q 012 ( 10 ^ 

4.8 

0.6 

0.6 

“ 011 in/ln/°F) 

2.2 

0.4 

16.4 

0022 00"® ln/tn/°F) 

2.2 

16.4 

0.4 


Note the 0-ply and 90-ply properties are complementary as expected. 
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(2) Th« reap«ctive thlcknenH ratios for this 8-ply APL are: 


Number of plies 




+45 

4 

4/8 

- 0.50 

0 

2 

2/8 

- 0.25 

90 

2 

2/8 

- 0.25 


(3) The corresponding Qc'* 

‘Jcxx " ''p+45^±451l * ''poQoil Vp9oQ90H 

- (0.50 X 6.1 + 0.25 X 18.7 * 0.25 x 2.0)xl0^ - 8.22x10** pal 


Qcyy ■ Vp^4 5^^14522 ''po‘Jo22 Vp9oQ9022 

- (0.50 X 6.1 ^ 0.25 x 2.0 + 0.25 x 18.7)xl0^ - 8.22x10^ psl 




^cyx " *^cxy ■ ''p±45‘^f4512 VpoQoU Vp9oQ9022 

- (0.5 X 4.8 -f 0.25 X 0.6 + 0.25 x 0.6)xl0^ - 2.70x10^ pal 


(4) The APL moduli and Poisson's ratio are: 


^cxx 

- Q 

^cxx 

Qcyy 

(<'•” - 8'.?2 , 

Jx 10^ - 

7.33x10^ 

psl 

Ej.yy 

■ O — 

Vcyy 

«?yx . 
Qcxx 

(-•” - M’f] 

1 X 10^ - 

7.33x10^ 

psl 

'^cxy 

■ 

■ 

iOfl 

8.22 


a 

0.328 



cyy 





''cyx 

- 

^cxx 

2.70 

8.22 


m 

0.328 



Check: vcyx - v^xy 

0.328 - 0.328 


> 0.328 - 0.328 O.K. 
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Th* calculations for F.(.yy« v^y^ and for the "check" were carried out for 
completencHK since it is obvious that for this loninate ^cyy " l^cxx 
''cyx ■ ''cxy 


(^) The TECs for this Al’L are calculated from equations (2.1) ind (2.2) but 
are generalized to include mure than two different ply combinations. The form 
of the equations using the summation sign are: 


“cxx 


±4S,0,90 


'»cxy‘^e21>®0ll * <‘^012 


''cxy‘^022)“02il 


“cyy 



+45,0,90 



'' cyx ^ 011^“011 ^^^022 " '' cyx '' 021^®023 

(2.9) 


where the sum is taken over 0 ■ ^^3, 0 and 90. Using the values calculated 
previously in equations (2.9) and cancelling the 10^ term with the 10"^ within 
the braces we have: 


7.33x10^ 


1 


|o.50[(6.1 - 0.328 X 4.8)(2.2) + (4.8 - 0.328 x 6.1)(2.2)] 

+ 0.25["(18.7 - 0.328 x 0.6)(0.4) + (0.6 - 0.328 x 2.0)(16.4)J 
+ O.25Q2.O - 0.328 X 0.6H16.4) + (0.6 - 0.328 x 18.7)(0.4)]| 

|o.50[9.96 + 6.1^ + 0.25[j.40 - 0.9^ + 0.25[29.6 - 2.2l]j 


°k:xx 


cyy 


7.33x10 
2.25x10-6 In/in/op 


^ r (o.SORa.S - 0.328 x 6.1)(2.2) + (6.1 - 0.328 x 4.8)(2.2)1 

7.33x10'’ I 

+ O.25Q0.6 - 0.328 X 18.7X0.4) + (2.0 - 0.328 x 0.6)(16.4)J 

+ 0.25[(0.6 - 0.328 x 2.0)(16.4) + (18.7 - 0.328 x 0.6X0.4)J | 

|o,50[6.16 + 9.9^ + 0.25[]-2.21 + 29.^ + 0.25[^0.92 + 7.4^| 


7.33x10' 


“cyy 


2.25x10‘6 ln/ln/«F 


I 


! 
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JIM vxpecled. The reader may f Ind it InMtructive to not* that: (1) the values 

wlrMn the various parentheses for ac»y and "cyy complementary and 

(2) the values for and <>cyy ' *hout equal to those for and 

“022 ® 

The reader will obtain valuable practice and insight by using the pro- 
cedure to calculate TECs of APL with ply configuration of his choice and a 
different composite system. 

3.0 TRANSFORMATION OF THERMAL EXPANSION COEFFICIENTS 

The thermal expansion coefficients (TECs) about any x'-y' coordinate 
axea of an orthotroplc angleplled laminate APL with material symmetry about 
the x-y coordinate axes are given by: 


“cx'x' 

" “cxxC°8^^ acyysin2^ 


"■cy'y' 

? 2 

■ a(.xxSln^^ + OgyyCOS ♦ 

(3.1) 

“cx'y' 

■ (<»cyy “ “cxx^“^° 2^ 



where the notation In equation (3.1) is as tollowa: Qcx'x'* “cy'y* 

“cx'y* TECs about the new coordinate system x'-y'; oexx ai'd ®cyy 

are the TECs about the x-y coordinate system and are calculated as described 
in Section 2; ^ Is the angle that the x' axis makes with the x axis. 

Note, Ogjj'y* is a shear-type thermal deformation which is present along any 
coordinate system x'-y' located at some angle ^ where 0 < ^ < 90 
since sine ^ 0 In this range. 

To perform the calculations using equations (3.1) we need the TECs Qexx* 
ucyy and the angle The TECs for the APL of Interest are either known or 

can be computed using the procedure and examples described in the previous 
section. The angle ^ Is known once the coordinate x'-y' axes has be *n 
selected. We can check our results using <^cx'x' “cy'y' " “exx ^cyy' 

Example 3.1. - Calculate the TECs of the [^0/o£| _ APL made from AS/E 

composite about an x'-y' coordinate axes where the x'-axls is located by 
^ ■ 15° from the x-axia. From example (2.1) oexx “ -0.68x106 in/in/°F and 
Ocyy ■ 13.0 x 10~6 in/ln/°F. Using these values in equations (3.1) we have: 

a^x'x' “ Ocxx°°®^^ OcyySin^^ 

- (-0.68 cos2 15° + 13.0 sln2 15°) x 10"^ in/in/°F 

- 0.24x10'^ ln/ln/°F 


(Cont 'd. ) 
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“cy’y* 


> a^yyco.24 


®cx’y’ 


®cx'y' 


(-0.68 Jln2 IS** 4 13.0 co- 15°) x ln/ln/°F 

12.08x10"* ln/ln/°F 

<“cyy ’ “cxx)*^*^® 

|[l3.0 - (-0.68)]sln2(15®)| 10"* ln/ln/°F 
6.84x10"* In/ln/op 


Aa can be seen, Che ahear-type TEC a^x'y' aubatantial. Reat raining thla 
APL along the x'-y* coordinate axea will induce conaiderablc In-plane ahear 
atreaaea. 


4. UkMINATE THERMAL STRAINS AND STRESSES 


Along the laminate m at erial axea . - The equal iona for calculating thermal 
atralna c^xy ^cyv alo''K the laiminaCc x-y coordinate axea are: 


*cxx ■ (T •* Tq)Ocxx 
tcyy ■ (T “ To)aj-yy 


(4.1) 


where T is the uae temperature and Tq ia the reference temperature (uaually 
room temperature); and c^-xy and acyy the TECa which are either known or 

can be determined as described previously. Use of equations (4.1) requires 
Chat Che TECs be independent of temperature within the range T - Tq* 

Example 4.1. - Calculate the thermal strains for the {+.30/0^5 APL made 
from AS/E composite (Example 2, 1) where T • 270° F and Tq ■ 70° F. The 
values for the TECs from Example 2.1 are: 

- -0.68x10"* in/ln/°F; a^yy - 13.0x10"* ln/ln/°F 
Substituting these values in equation (4.1). we obtain: 


^exx " <T - To)acxx “ (270 - 70) (-0. 68x10"*) - -136x10"* In/ln 
Ccyy - (T - To)Ocyy - (270 - 70) (13.0x10"*) - 2600x10"* In/in 

The equations to calculate the corresponding thermal stresses, assuming Chat 
che laminate is completely restrained from thermal expansion, are: 

°cxx “ "^T'(Qcxx“cxx Qcxy®cyy) 

(4.2) 

°cyy “ “AT(Qcyxacxx + Qcyyocyy) 
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whrr* AT • (T • T^)); th« Q^'m and «r« d»t«.'rain«d aa deacrlbad 

prcvioualy . 


E xaapl ^.2. - Calculate tlia raatrained ttienaal atrraMea In Exaapla 4.1. 
Kaf erring to Rxamplaa 2.1 and 4.1, we have: 


- 15.0x10* pal 
Qcyy ■ 2.45xlO*» pal 
^cvx ■ ‘Jcxy • 2.15 x10<> pal 
Ualng theae vaiuea In equations (4.2) 


®cxx " in/ln/OF 

Ocyy ■ 13.0xlO-<> in/in/OF 
AT - (270 - 70) • 200‘’ F 
we obtain (cancelling 10^ with 10~^) t 


cxx 


’cyy 


*T(Qcxx®cxx * % xy \ yy ^ 

- 200[l5.0 X (-0.h8) ♦ 2.15 x 13.()] 

- 200(-l0.2 ♦ 28.0) 

- 3550 pal - -3.6 kal 

- AT(Q a 4- Q a ) 

' *cyx cxx ^cyy^cyy' 

- 200[2.15 X (-0.6C) + 2.45 x 13. oj 

- 200(-1.46 + 31.8) 

- 6078 pal ■ -6. 1 kal 


Two points are worth noting In connection with the above values of these 
thermal stresses: 

1. The thermal stresses Ocxx Ocyy ®re relatively small (4 and 28 

percent, respectively) compared to the corresponding compressive failure 
stresses of the l.imlnate ■ 83 ksl and S^yy^ ■ 22 kal, based on first 

ply failure (ref. 5). 

2. The magnitude of these thermal stresses may be sufficiently high to 
cause panel buckling. For example, a 20 In. x 10 In. x 0.04 In. panel from 
this 8-ply APL has buckling stresses of about Oexx * ^^0 psl and 

0(.yy * 920 psl (calculated using the equation in ref. 6) or approximately 15 
percent of the restrained thermal stresses which are relatively low. A panel 
with this geometry will buckle at an Increase In temperature of about 15 per- 
cent of 200^' F or 30° F. The Important conclusion from the discussion in this 
example Is that thermal stresses need be considered carefully by the designer/ 
analyst In situations where restraints may be present. 


5.0 PLY THERMAL STRAINS AND STRESSES 

It Is Instructive to describe the ply thermal strains and stresses in the 
plies of an APL by breaking them down into four "commonly thought-of" types. 
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These types ere: 

1. Restrained - APL is restrained froe thermal ex*>unsiun. 

2< Free - APL is free to undergo thermal expansion. 

3. Residual • APL laminate Is cooled from cure temperature to 

use temperature and frequently to room temperature. 

4. Combined - Combinations of free and residual. 


The ply thermal strains and stresses to be described are those slong the ply 
material axes 1, 2, 3« figure 1. For convenience, the strains along the 
1-dlrectlon (fiber direction) are defined by C(ix stresses otxil 

those along the 2-directlon (transverse to the fiber direction) are defined 
by C|22 *’122* those In the 1-2 plane (Intralaminar shear) are de- 
fined by cti2 °tl2* 

Restrained aPL . - The ply thermal strains for this case are given by 

«112 ■ ^122 • till ■ 0 (5.1) 


The corresponding ply stresses are given by 

• ■ *T(Qtuaixi + Qll2“t.?2) ^ 

®122 ■ “ ^T(Q£21«U1 * ‘^I22“t22> f 

°tl2 " ® > 


(5.2) 


where AT equals the use temperature minus the reference temperature. Where 
the Qt's are the reduced ply stiffnesses and the oi's are ply thermal 
expansion coefficients (ply TECs). The ply reduced stiffness Q{,'b and TECs 
can be estimated from figures 3 to 18 at b ■ 0*’. Equations (5.2) show that 
the ply material axes thermal stresses In an APL restrained from thermal ex- 
pansion depend only on ply properties. Also, there Is no intralaminar shear 
stress for this case. 


E xample 5.1. - Calculate the ply thermal stress in the piles of the 
[^30/O^g APL, made from AS/E composite, where the temperaturi Is increased 
from 70® F (room temperature) to 270° F and the APL Is restrained from thermal 
expansion. The numerical values we need for these calculations are deterroinr t 
as follows: Figures 7 and 8 at 8 ■ 0° yield 

Qtll “ 18-7x10** psi, Q^22 • 2.0x10^ psi 


Qj 21 " *^112 ■ 0*8xl0** psi 


ojii • 0.4 x10"8 in/in/°F; 0£22 * 16.4x10-8 in/ln/°F 
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«nd AT ■ 270° F - 70“ F • 200° F coniilatent with the prevloua definition. 
UelnR these numericsl values in equations (5.2), we calculate 

‘’ill “ " ^^^‘^tll‘*tll ‘^I12”l22^ 

- - 200(18.7 X 0.4 0.6 x 16.4) - - 3464 pal • - 3.5 ksi 

°122 ■ ^'‘'^‘^t21°tll '^I22“l22^ 

- • 200(0.6 X 0.4 2.0 x 16.4) - - 6608 psi - - 6.6 ksi 

Free APL. - The ply thermal strains for this case are given by: 

Cm - ATta^^^cos^e + u^yysln^e - 0^^) 

^122 " AT(aj,^x*^°^® Ocyy^®*^® " °122^ (5.3) 

^112 * 6T(agyy - aj.jjjj)sin 9 


The corresponding ply stresses are given by: 


°lll • ^T[Qiii(aj.,^^co820 acyy‘*^"^®-‘'tll> +‘Jll2<“cxx»i"^® + “cyy‘'°®^®'“l22)] 

°122 “ ^t[Qi2i(«cxx°°‘*^® “cyy®^"^®'“Ul)+Ql22(“rxx®^"^® + “cyy‘=°®^®-“l22)] 

0112 “ ATi;Qj 33 (acyy - a^^^lsin^oj 


(5.4) 

Also, when the ply thermal strains are calculated using equations (5.2), the 
corrr :>ponding ply stresses are given by: 

°lll “ ‘^lll'^lll ^112*^122 

°122 “ ‘^I21‘=lll ^122^122 ^^-5) 

°112 “ ^133^133 

where the strains (e^) are calculated from equation (5.3). 

_ E:^mple 5.2. - Calculate the ply therro-il strains in the plies of the 
[± 30 / 02 ] s APL, made f rom AS/E composite, where the temperature is increased 

from 70<^ F (room temperature) to 270° F and the APL is not restrained from 
thermal expansion. We will calculate the ply thermal strains in the ±30° 
plies, in the -30° plies and in the 0°-plles by using equations (5.3). The 
numerical values we need are: 


lA 


AT ■ 200'* F 

0 - lO**, -30°. 0“ 

®cxx " ln/in/”F ^ 

a^yy - 13.0x10"** ln/ln/°F 
“til ■ 0‘^xlO"* ln/ln/“F 
0122 • 10-AxlO"** in/ln/^F ^ 


(Theite values arc taken from Example 2.1) 


3 0°-ply . - SubatltuCing these numerical values and 0 • 30° In equa- 
tions (5.3) we calculate: 

"ill “ ^T<“cxx"°“^« ^ “cyy«^"^® " “tll> 

- 200(-0.68 X cos 230° 13.0 x sln230° • 0.4) x 10"^ in/tn 


• 468x10”** in/in or ■ 0.05Z 

"122 “ AT(aj.jjjjSln^0 + aj.yycos20 _ 

- 200C-0.68 X 8in230° + 13.0 x co8230° - 16.4) x 10"* in/in 


- - 1364x10"** In/ln • -0.14X 

"tl2 • *T(a^yy - acxx>“^"^® 

- 200[l3.0 - (-0.68)] X 10'^ sin 60° 

• 2369 X lO"^ in/in • 0.24Z 


-30°-ply . - The thermal strains in the -3Qo ply are the same as those in 
the -♦■30° plies for ejn and £^22 opposite sign for Cfi.2* The reader 

can readily verify that this is the case by inspection of the appropriate 
equations. 


0°-ply . - Substituting the above numerical values and 0 
tions (5.3) we calculate: 


0 


o 


in equa- 
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f|lU “ '*cyy**^'* ® " ‘*11 1 ^ 

- 200(-0.68 X coB^0“ * 13.0 X Bln^O** - 0.4) x 10“^ in/ In 

- -216x10"** In/ln - - 0.C2X 

2 2 

*122 ■ AT(iij,xx**^“ ® ■*■ aryv**'** ® " ®t22^ 

- 200(-0.68 Bln^O® + 13.0 r con^O° - 16.4) x lO’^ 

- -680x10"*’ In/ in " - 0.07Z 

cti2 " AT(c*j.yy “ apxx^**^*' 

- 200|l3.0 - (-0.68)) X lO'** sin 2(0°) 

■ 0.0 


Tin* ri'ader will t lud It InBtrui'tlve t«> oompare the corresponding thermal 
strains In the 30*’ and O" piles. Both normal strains C£H and €£22 the 
30*’ ply are about twice those In the 0‘* piles wlille the shear strain E£j^ 2 
has the largest magnitude In tlie 30*’ piles and Is "zero" In the 0*’ piles. 

E xamp le - Calrulate the corresponding thermal stresses In the piles 

In the API. of Example 5.2. To calculate tlie corresponding thermal stresses 
we can use either equations (5.4) or (5.5) since we already calculated the 
strains In Example 5.2. We will use equations (5.5) for convenience. In 
order to use equations (5.5), we need numerical values for tlie reduced ply 
stiffnesses ()£ and corresponding thermal strain values l£ from Example 5.2. 
We tabulate these numerical values for convenience. 


Reduced Ply Stiffnesses In 10® psl 
(fig. 7 at 0 ■ O*’ or fri>m 
Example 5.1) 

Q£11 ■ 18.7 

Q£22 “ 

*U21 “ *Ul2 “ 0-60 

Qt33 “ 0.5b 


Thermal Strains In 10”® In/ In from 
Example 5.1 for 9 ■ 



30*’ 

-30*’ 

0*’ 


468 

468 

-216 

*t22 

-1364 

-1364 

-680 

*U2 

2369 

-2 369 

0 


Using corresponding values in equacions (5.5) we have (cancelling 10 with 10 
for convenience) 

-iPly • Olll • Qtlicill Qtl2Ctl2 


- 18.9 X A68 + 0.6 x (-1364) - 8367 psi 


°122 “ ‘^121‘^m ^ ^122*^122 


- 0.6 X 468 2.0 X (-1364) - -2447 psi 


0^-ply ; 


®tl2 " ^633^133 " 0.56 x 2369 ■ 1327 psi 
* ®tll ®t22 same as for the •♦•30^ plys 

an2 has opposite sign or 0ji2 * " 1327 psi 

°ui “ ^ui‘=ui * 


- 18.9 X (-216) + 0.60 X (-680) - -4490 pal 
°222 “ ^ 121^111 ^ 122^122 


• 0.60 X (-216) + 2.0 X (-680) - -1490 psi 


Oj 22 " ^£33*"£33 " 0,56 x (0) “ 0 


The interesting point to be noted from the numerical values of the ply 
thermal stresses is that the transverse ply stresses 0|^,22 compressive 

and are only about 22 percent of the compressive ply strength (S ^220 
about 20 ksi at 270° F). This implies that thermal fatigue in the range 
70° < T ^ 270° is not anticipated to cause progressive degradation to the 
ply transverse properties in the A**! considered. 


6.0 PLY LAMINATION RESIDUAL STRAINS AND STRESSES 


Ply lamination residual strains and 
thermal strains and stresses. However, 
posite laminates they need be treated in 
and stresses arise from the fabrication 
"The difference between the cure and use 
They are always present in the plies of 
magnitude need be determined to accurate 
rity of APLs in structural applications. 


stresses are a special case of 
because of th“ir Importance in com- 
a separate section. These strains 
procedure of the composite laminates: 
temperatures" as mentioned previously 
APLs (free or restrained) and their 
ly assess the thermomechanical Integ- 
The laminations residual strains 
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can be calculated usIdk t'(|uatlunH (S.3) and tb<’ laminatiuii realdual Htrt*HMi‘« 
can be calculated ualng equatlona (5.4) or equatlona (5.5) when the lamlna- 
tlona residual straina are either known or they have been calculated pre- 
viously. Usually, the lamination residual strains and stresses are calcula- 
ted at room temperature. The temperature difference AT for these cases is 
known from the specified cure conditions. For example, the AT for APLs 
made from structural epoxies is about -300° F. The minus sign is determined 
from the definition of AT given In section 5.0: 

AT • RiH)m lempcraturc* - Cure Temperature 
(about 70°F) (about 170*’ K for 

structural epoxies) 


Calculations of lamination residual strains and stresses will be Illus- 
trated using the Al’l.s considered In Kxample 5.2. 


Example fi.l. - Calculate the lamination residual strains and stresses In 
the piles of the [+30/0^18 API, made from AS/K composite. First, we calculate 
the lamination residual strains using equations (5.3). The numerical values 
we need to use in these equations are obtained from Examples 5.2 and 3.3 and 
are summarized below for convenience. 


a - -O.bdxlO"^ ln/ln/°F 
exx 

■ 1 1.0x10“^ ln/ln/°F 
®tll “ 0.4x10”^ ln/ln/°F 


■ 18.7x10^ psi 
'^22 “ 2.0x10^ psl 
*^121 * ^tl2 " 0.60x10^ psl 


“ lb. 4x10*^ ln/ln/‘’F “ 0.5bxl0^ psl 

0 - 30°, -30°, 0°; AT - (70° F - 370° F) - - 300° F 

I'slng appropriate numerical values In equations (5.3) we calculate the lamin- 
ation residual strain: 


+30° PLY: 
0° n.Y: 


AT(a^,^x‘'°«^® + u^yyHin^e - aj,j) 

-300f(-0.68) X 0.75 + 13.0 x 0.25 - 0.^ x lO"^’ - -702x10"^ in/ln 
-100[(-0.b8)(l 0) + 1 1.0 X 0.0- 0.4*] x 10"^ -• 324x10"^ 


In/ in 
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ei22 

Pl.Y; 

0*’ PLY: 

*^112 

±30‘’ PLY; 

0° PLY; 


2 

AKog^^sln^e + OcyyCOB 0 - ai22^ 

-l(K)[(-0.68) X 0.25 n.O X 0.75 - 16. a| x 10"* • 2046x10"* In/ln 
- 300 Q-O. 68 ) X 0.0+ 13.0 x 1.0 - 16. 4j X 10"* - 1020x10"* In/ In 

AT(Oj.yy - 20 

- 3 O 0 Q 3 .O- (-0.68)^ (± Bin 60”) X 10"* - + 3554xl0“^ tn/ln 
-300[l3.0- (-0.680(0.0) X 10"* - 0 


UsinR appropriate ()£ values and C£ values In equations (5.5) we calculate 
the ply lamination residual stresses (cancell ing 10* with 10“*): 


t30” PLY: 
0 ” PLY: 


°tll “ ^tll‘-£ll ^ 112^122 

• 18.7 X (-702) + 0.60 x 2046 - - 11 900 psi • - 11.9 ksi 
- 18.7 X 324 + 0.60 x 1020 - 6671 psi • 6.7 ksi 


®£22 “ ‘^£ 21‘^£11 ‘^£ 22‘^£22 

+30” PLY: - 0.60 x (-702) + 2.0 x 2046 - 3671 psi - 3.7 ksi 

0” PLY: - 0.60 x 324 + 2.0 x 1020 - 2234 psi ■ 2.2 ksi 


*^£12 “ ^£ 33*^£33 

±30” PLY: - 0.56 X (^ 4104) - + 2298 psi - + 2.3 ksi 

0° PLY: - 0.56 X 0.0 - 0.0 


The transverse lamination residual stresses in the ±30 plies are about 50 ner- 
cent of the corresponding ply strength (6 to 8 ksi). 




7.0 CONCLUDING RKMARKS 

A convenient procedure !■ described to determine the thermal behavior (ther- 
mal expansion coefficients, thermal and residual stresses) of angleplied fiber 
composites. Tlte procedure consists of equations and appropriate graphs for var- 
ious (^6) ply combinations . 'Dtese graphs consist of reduced stiffness and ther- 
mal expansion coefficients for frequently used composites and hybrids as func- 
tions of X.6 in order to simplify and expedite the use of the equations. The 
procedure is applicable to all types of balanced, symmetric fiber composites in- 
cluding interply and intraply hybrids. The versatility and generality of the 
procedure is Illustrated using several step-by-step numerical examples. The step- 
by-step numerical examples are set up so that the calculations can be made using a 
pocket calculator. Sooie of the numerical examples ware selected to illustrate 
significant implications of composite thermal behavior in design applications. 


APL 

AS/E 

B/E 

E 

c 

HMS/E 

K/E 

Qc 


Qe 


s 


I 


S-C./E 


T 

T 

o 

AT 

TEC 


8.0 SYMB01.S 

angleplied laminate 

AS-graphlte-f Iber/epoxy-mat r lx compoHlte 
boron-f Ibcr /epoxy-matrix compoalte 

laminate modulus - subscripts x,y denote structural axes directions 
ply modulus - subscripts 1,2 denote ply material axes direrr one 
high modulus graphite-f iber/epoxy-matrlx composite 
Kevlar-f Iber/epoxy matrix composite 

reduced l.imlnate stiffness - subscripts x.y denote structural axes 
directions 

reduced ply stiffness - subscripts x.y denote ply nuiterial axes 
directions 

reduced stiffness for tU symmetric laminate - subscripts 1,2 
denote material axes directions 

ply strength - subscripts 1,2 denote ply material axes directions; 

- subscripts T, C, S denote type 

S-glass-1 Iber/epoxv-raatr lx composlt e 

temperature 

reference temperature 

temperature difference between use and reference temperature 
thermal expansion coefficient 
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Vp ply thickntiiN ratio - aubacripta 0, 0, 90 dt*noCe ply dca Ignat Icm 

to which thi' ratio appllva 

x,y.x atructural axea coordinate dlrcctlona 

1,2,3 aatcrlal axea coordinate dlrcctlona • 1 taken along the f iber direct Ion 


laminate configuration dealgnatlon - numbera In the blanka denote 
ply aracklng aequence and orientation - aubacrlpt S denotea 
aymmetry about ply In laat blank a^ace 

Oc laminate TEC - aubacrlpt a x,y denote laminate structural axea 

direct Iona 


C-/./-1 


“l 

“0 




ft 

6 

V 

c 


ply TEC - aubacrlpts 1,2 denote ply material axea dlrectlona 

laminate TEC - aubacrlpta 1,2 denote material axea dlrectlona 

laminate atraln - aubacrlpta x,y denote atructural axea dlrectlona 

ply atraln - aubacrlpta 1,2 denote material axea dlrectlona 

ply orientation angle meaaured from the x-lamlnate atructural axea 
to the 1-ply material axes and taken positive 

laminate Poisson's ratio - subscripts x,y denote atructural axes 
directions 


ply Poisson's ratio - subscripts 1,2 denote ply material axes 
directions 


♦ 


laminate stress - subscripts x,y denote structural axes directions 

ply stress - subscripts 1,2 denote material axes directions 

laminate coordinate axes x', y', z' orientation other than the 
structural axes x, y, z measured from the x axis to the 
x'-axls and taken positive. 


Conversion factors: 


”C - |(°F -32°) 

A°C - I A °¥ 

cm/cm/“C • "I in/in/”K 


MPa 


b.89 ksi 
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